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1 Introduction 



1.1 Motivation 



This 



paper establishes a stochastic integration theory for volatility modulated Volterra (VA4CV) 
processes which are defined as 



X(t) = f g(t,s)a(s)dL(s), 
Jo 



(1) 



where g is a deterministic function, a is a predictable stochastic process and L is a Levy 
process. The precise regularity assumptions needed to ensure that the integral in (pQ) is well- 
defined are given in the following Subsection. 

O ur rese arch is closely related to the influential work bv lAlos et al. I (|200lh and iDecreusefondl 



(|2002l . 120051 ). who developed a stochastic calculus for Gaussian processes. Our new contribu- 
tions include in particular a stochastic integration theory with respect to Volterra processes 
which include a stochastic volati lity /interm i ttency component, which is denoted by a in ([1]). 
As such we extend the results bv lAlos et al.l ( 200ll ) in the Gaussian case to allow for volatility 
modulation. Moreover, we also develop a suitable stochastic integration theory for the case 
that L is a pure jump Levy process. Also these new results hold in the presence of stochastic 
volatility. Furthermore, ou r proposed integral is of a slightly different form than the one sug- 
gested in lAlos et al.1 (j200lh . In addition to developing the new integration type, we establish 
the main properties of our new integral, and we focus on a variety of interesting examples. 
E.g. we study integration of deterministic kernel functions with respect to Volterra processes 
and we investigate the case when the integrand is itself a Volterra process, which leads to a 
representation result in terms of the second chaos expansion. 

Due to the flexibility the class of VM.CV offers, it can be applied to modelling a wide 
range of phenomena. E.g. VA4CV proc esses have recently been found to be suitable for 
e.g. modelling turbulence in physics, see Barndorff-Nielsen fc Schmiegel ( 20091 ). Also, they 
can be used to model financia l data, such as commodity spot prices, as recently shown in 
Barndorff-Nielsen et al.l (|2010l ). In many applications, it is important to have a stochastic 
integration theory at hands, e.g. in financial applications one typically models the value of a 
portfolio by a stochastic integral with respect to the price process. This paper aims at estab- 
lishing the suitable integration theory which can then be used in various such applications. 

The outline for the remaining part of the paper is as follows. First we review the technical 
assumptions needed to define integrals of the type (pQ). Next, since Malliavin calculus both 
for Brownian motion and for pure jump Levy processes is the key tool for establishing our 
integration theory w.r.t. VM.CV processes, we give a short review of main results from that 
calculus in Section [2j Section [3] introduces the class of volatility modulated Volterra processes 
in more detail and describes relevant examples. The main contributions of the paper are 
contained in Section U where we establish the integration concept with respect to volatility 
modulated Volterra processes - both in the case when the process is driven by a Brownian 
motion and also in the case of a pure jump Levy process. Next we derive the key properties 
of the new integral in Section [5] and we focus on various applications of the new concept. 
Finally, Section [6] concludes and gives an outlook on future research. 
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1.2 Preliminaries 



Let (fi, J 7 , P) be a complete probability space equipped with a nitration # = (J r t)j>o satisfying 
the usual conditions of right-continuity and completeness. 

We aim to develop a stochastic integration concept with respect to volatility modulated 
Levy-driven (VA4CV) processes. In order to set the scene we first review the suitable in- 
tegration concept for proper integrals with respect to a Levy process L, i.e. integrals of the 
form 

<Kt)dL(t), (2) 



f 

Jo 



where <f> = (<f>(t))teR is a predictable integrand. Clearly, by replacing 4>(t) by (j>(t)I{ a <t<b} for 
0<a<6<oowe are able to treat integrals over th e intervals ( a, b] w it hin this framework 
as we ll. Relevant references in this context i n clude iKallenberd (119751) . iRajput fc Rosihski 

(ll989h.lKallenberg| (Il992l] jDecreusefondl (l2002h . IJacod & Shirvaev (120031. Theorem III.6.30 ), 

Chernv k Shirvaev! (|2005 h . IPi Nunno et all (|2005l l. |Protterl (|2005l ). lNualartl (j2006Eh , basse- O'Connor et al 
(l2(U2h . 

Let us recall some key properties of Levy processes which we will apply throughout the 
paper. 

Definition 1. A stochastic process L = (L(t))t>o is a Levy process if it satisfies the following 
five conditions: (1) £(0) = a.s.; (2) L has independent increments in the sense that for 
any choice of n > 2 and t\ < ■ ■ ■ < t n the increments Lfo) — L(ti), . . . , L(t n ) — L(i n _i) are 
independent; (3) L has stationary increments in the sense that for all s < t, L(t) — L(s) has 
the same law as L(t — s); (4) L is continuous in probability; (5) L has cadlag paths. 

The Levy process L has a Levy-Khintchine representation, which is given by 



E(exp(iuL(i))) = exp(tip(u)), ip(u) = iju 



c 2 u 2 



+ 



l-iuT{z))i{dz), (3) 



where r is a truncation function which is assumed to behave like the identity function around 
the origin. More precisely, a standard choice would be t(z) = zl^^iy. Further, 7 G R, 
c 2 > and t denotes the Levy measure of L satisfying ^({0}) = and L(l A \z\ 2 ) £(dz) < 00. 
We write (7, c 2 ,£) for the characteristic triplet of L. Recall that we have the following Levy- 
Ito-representation for < s < t: 



L(t) - L(s) = (t- s) 7 + c(B(t) - B(s)) + 



r(z)N(dz,du) 



(4) 



+ 



/'/< 



z — t(z)) N(dz, du), 



where N denotes the corresponding Poisson random measure and N(dz, du) = N(dz, du) — 
£(dz)du denotes the compensated Poisson random measure. Moreover, B is a standard Brow- 
nian motion which is independent of N . 



1.2.1 Stochastic integration with respect to Levy processes 

Next we review briefly how the integral in (ITT) is defined. For de t ermin istic integrands one can 
apply the integration concept developed by Raiput Sz Rosihski ( 19891 ) for Levy bases as inte- 
grators, using the fact that Levy processes can be directly linked to homogeneous Levy bases, 
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which are def ined as homogeneou s , infin i tely d i visibl e, independently scattered random mea - 
sures, see e.g. iRajput &: Rosihski (|l989h . lSatol (120041 ) . [Ba7ndorff-Nielsen fc Schmieeel (120071 ). 

The integration theory can be extended to allow for stochastic integrands as follows. 
Recall that £ = (•? r t)t6K denotes a filtration satisfying the usual conditions. Let V denote the 
predictable cr-field on R x Q, i.e. 

V = a ((s, t) x A : -oo <s<i<oo, J 4e 7" s ) 

which is the smallest sigma-field generated by the left-continuous and adapted processes. 

Assume that L = (Lt)t>o is a Levy process with respect to 5 with characteristic triplet 
(7,0-V). That is, 



for all < s < t, L(t) — L(s) is independent of F s 



(5) 



As shown in iBasse-O'Connor et al.1 (j2012l ). see also IChernv fc Shirvaevl (120051 ). the integral 
Jo°° 4>{ s )dL(s) exists if and only if <f> is predictable and the following three conditions are 
satisfied a.s.: 



a 2 (j) 2 (s) ds < 00, 



[1 A {(f>(s)z) 2 ) £{dz)ds < 00, 



*)7+ / [r{4>{s)z)-cp{s)T{z)\^dz) 



(6) 



ds < 00. 



1.2.2 Important examples 

We will now focus on some important examples, where the integrability conditions stated in 
© simplify considerably. 

Example 1. If in addition to ©, for all s < t, L{t) — L(s) has zero mean, is square-integrable, 
and the predictable process <p satisfies E( (p(s) 2 ds) < 00, then ^ exists, is square- 



integrable and Jq 4> s dL s is a square-integrable martingale up to infinity, see lBasse-O'Connor et al 
(|2012l ) for more details. 

Example 2. Suppose L is a Brownian motion with drift. Then the characteristic triplet 
is given by (7,<r 2 ,0) and condition ([5]) is satisfied. Then the integrability conditions for a 
predictable process 4> simplify to 



(s) 2 ds < 00, 



/ \4>(s)\ ds < 00 a.s. 
Jo 



Example 3. See IChernv fc Shirvaevl (120051 . Corollary 3.6) for a detailed treatment of the 
case when (L t )t>o is an a-stable Levy process. 



2 A brief background on Malliavin Calculus 



We briefly review the key results from Malliavin Calculus needed to establish our new stochas- 
tic integrat i on the ory. I n the following, our pr esentation is extracted from the monographs 
by iNualartl (|2006^ and bi Nunno stall (|2009h . where we discuss the Wiener and the pure- 
jump case separately. 
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2.1 The Wiener case 



As before, (fi, T, P) denotes a complete probability space and we let (T,B,/i) be a measure 
space with fj, a cr-finite measure without atoms. The space of square-integrable functions on 
(T,B,fj) is denoted by L 2 (T). In our applications, T is typically the finite time interval 
T = [0, T], for some given finite time horizon T < oo, equipped with the Lebesgue measure 
fM on the Borel sets B. 

We denote by W the L 2 (f2) Gaussian measure on (T,B). That is, for every A G B with 
fJ>(A) < oo, W(-A) is a centered normally distributed random variable with variance n{A). 
Furthermore, if A, B G B with /j,(A) < oo, fJ,(B) < oo, are disjoint, then W(.A) and W(B) are 
independent. Assume that T is the cr-field generated by W(A) for A £ B with //(A) < oo. 

To define multiple stochastic integrals with respect to W, consider elementary functions 
/ € L 2 {T n ,B n ,n n ) given as 

m 

f(ti,...,t n )= a>i 1 —i n 1-Ai 1 x-xAi n (tl}---,tn)j 
h,— >*n=l 

for Ai j G ^ such that /x(A^) < oo and with the property that a^...^ = if any two of the 
indices ii, . . . ,i n are equal. Then, we define 

m 

W)= E a il ... in W(A il )---W(A in ). 

il,— ,i„=l 

The set of such elementary functions is dense in L 2 (T n ) and the operator I n can be extended 
to a linear and continuous operator from L 2 (T n ) to L 2 (P). It has the properties that I n (f) = 
I n (/) if / is the symmetrization of /, and 

E[J m (/)J n ( 5 )] = { °' ~ 

I n\{f,g) n , m = n, 

where we have denoted the inner product in L 2 (T n ) by (-,-) n - The notation 

Uf)= [ f(h,...,t n )W(dt 1 )---w(dt n ), 

will be frequently used. 

The famous Wiener-Ito chaos decomposition (jNualartl ()2006a . Theorem 1.1.2)) holds for 
the space L 2 (P): 

Theorem 1. If X G L 2 (P), then there exists a unique sequence {f n }^Lo of symmetric func- 
tions where f n G L 2 (T n ) such that 

oo 

X = Y J In{fn)- 

71=0 

Moreover, 

oo 

\\X\\ 2 :=E[X 2 ] = Y,nU 



2 



where \ ■ \ n is the norm in L (T n ). 



n=0 

2, 
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The Malliavin derivative of a random variable X € L 2 (P) can be characterized as an 
operation acting on the chaos. The domain of the deriva tive can also be characterized by the 
chaos functions. The following discussion is taken from iNualartl (|2006a . Proposition 1.2.1). 
Suppose X € L 2 (P) with chaos expansion 



x = Y j i n {f n ) 



n=0 



is satisfying 



^2nn\\f n \l < oo. 

n=l 

Then we say that X € D 1 ' 2 , the domain of the Malliavin derivative Dt, and it holds that 

oo 

D t X = Y j m n -i{fn{;t)). (7) 



n=l 



The Malliavin derivative is thus the operation where we fix the last coordinate of the chaos 
function f n , and moves it to a chaos of order one n — 1. It corresponds to the annihilation 
operator in quantum physics. A direct computation shows that 



E 



(D t X) 2 n(dt) 



T 



^ran!|/ n ,| 2 . 



(8) 



n=l 



In particular, this proves that Dt is a linear operator from B 1 ' 2 into L 2 (Q, x T). 

The Skorohod integral is introduced as the adjoint operator of Dt in the following way: 
denote by Dom(5) the set of all processes Y in L 2 (Q, x T) where 



IE 



J Y(t)D t X n{dt) 



< c\\X\ 



0) 



for all X € B 1 ' 2 , and c > being a constant depending on Y . Then, for Y £ Dom(5), we 
define the Skorohod integral of Y, denoted by 



Sr 



Y(t)5W(t) , 



as the unique element of L 2 (P) characterised by 



E 



X J Y(t)6W(t) 



E 



T 



Y(t)D t X n(dt) 



(10) 



(11) 



for all X G D 1,2 . 

One may view the Skorohod integral as a creation operator on chaos, as the following 
shows. It holds by the Wiener-Ito chaos expansion that for any Y E L 2 (ft x T) there exists 
a sequence of symmetric functions f n (-,t) in L 2 {T n ) such that 



Y(t) = J2 Wn(; t)) 



n=0 
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In this representation, t is only considered as a parameter of the function f n (-,t) on T n . 
Letting f n G L 2 (T n+1 ) be the symmetrization of f n (-,t), it holds that 

r. OO 

/ Y(t)SW(t) = ^In + l(fn). (12) 

Moreover, the space of Skorohod integrable processes Dom(<5) is characterized by processes 
satisfying 



ln+1 < 0° 



5>+i)!|/ n 

n=0 

where this sum is equal to the variance of the Skorohod integral. 

The Skorohod integral is a generalization of Ito integration in the case of T = [0 , T] and 



[i being the Lebesgue measure on the Borel cr-algebra B. We have from iNualartl (|2006£ 
Proposition 1.3.4) that 

Y(t)SW(t) = [ Y(t)dW{t), 



'0 Jo 

whenever Y is Ito integrable. The essential extension offered by the Skorohod integral is that 
one can define the stochastic integral of Y with respect to Brownian motion without restricting 
to integrands Y being adapted to the filtration generated by the Brownian motion. Sometimes 
one refers to the Skorohod integral as an anticipative stochastic integral for this reason. 

Skorohod integration and Malliavin differentiation have several interesting properties, and 
we present the ones which are relevant for our analysis of stochastic integration. First, the 
following "fundamental theorem of calculus" holds: 

Proposition 1. Suppose that Y G Dom{5) flD 1 - 2 and j T Y(t)5W(t) G IS) 1 ' 2 . If the process 
s i— s> D t Y(s) for s G T is Skorohod integrable for almost every t G T , then 



A ^Y(s)5W(s^j = Y(t) + J^D t Y{s)5W{s) 



(13) 



for almost every t G T . 



The next result, which is an integration by parts formula, will be the key to define the 
stochastic integrals we are interested in this paper: 

Proposition 2 (Integration by parts formula). Suppose that Y G Dom(5) and X G D 1 ' 2 . // 
XY G Dom(5), then 

j XY(t)5W{t) = X i Y(t)5W(t)- j Y(t)D t X n(dt). (14) 



This result is proved and discussed in INualartl (12006 a . page 40) 



2.2 The pure-jump Levy process case 



Our ex position of the Malliavin Calculus for Levy processes is based on the monograph lDi Nunno et al 



(|2009h . Here, the space (7~, B, fi) is explicitly chosen to be T = [0, T] with fi being the 



Lebesgue measure on the Borel subsets B of [0,T]. We work under precisely the same as- 
sumption in the following review. 
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Suppose that is a pure-jump Levy process which is square integrable and with mean 
zero. Let N(dz,dt) denote the Poisson random measure associated to on [0,T] x M.q, 
with M.q := M\{0}. Moreover, we denote by £(dz) the Levy measure of L^ p \ In terms of its 
Levy-Ito representation, we have 

= [ [ zN{dz,ds), (15) 



where N(dt, dz) is the compensated Poisson random measure. For < t < T, we denote by 
J~t the cr-algebra generated by L^(s) for s < t, augmented by the sets of P-zero probability. 
Assume that T = Ft- 

Consider now the space L 2 (([0,T] x R ) n ,(dt x £{dz)) n ), that is, measurable real-valued 
functions / on ([0,T] x R ) n such that 

\f\ni ■= / f 2 (h, *uh, z 2 , ■ ■ ■ ,t n , z n ) dti £{dzi) dt 2 £(dz 2 ) ■ ■ ■ dt n £(dz n ) < oo . 

J{[0,T]xR ) n 

Furthermore, we denote by / the symmetrization of / in the pairs of n variables (ti, z\), . . . , (t n , z n ). 
Any function which is symmetric in these pairs is called a symmetric function from now on. 

The re-fold iterated stochastic integral of a symmetric function / G L 2 (([0, T] x Mo) n ) with 
respect to the compensated Poisson random measure is now defined to be 

I n (f):=n\[ [ ••• f [ f(t 1 ,z l ,...,t n ,z n )N(dz u dt 1 )---N(dz n ,dt n ). (16) 



JO JRo JO JR 

We have the Wiener-Ito chaos expansion for all square integrable random variables: 

Theorem 2. If X £ L 2 (P) is Ft -measurable, then there exists a unique sequence of sym- 
metric functions {fn}^Lo where f n S L 2 (([0,T] x Mo) n ) such that 

oo 

X = Y J Ufn). 

n=0 

Moreover, 



1*11! = E n! l^H 



n=0 

As in the Wiener case, we introduce a Malliavin derivative as an annihiliation operator on 
chaos. Since the kernel functions f n now are functions of pairs (t, z), it is natural to introduce 
a derivative operator D t , z parametrized over the pair (t, z) rather than t as was the case 
for Brownian motion. We say that X £ ID 1 ' 2 , the domain of the Malliavin derivative Dt >z , 
whenever 

oo 

^2 n n ] -\fn\ 2 n/ < OO. 

n=l 

In this case, we define 

oo 

D t>z X = Y j nI n -i{fn{;t,z)). (17) 

n=l 
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Again, as in the Wiener case, we can introduce a Skorohod integral as the adjoint of the 
Malliavin derivative. We say that the space-time random field Y € L 2 ([0,T] x Ro x Q) is 
Skorohod integrable (that is, belongs to Dom(£)) if 

cT 



|E 



Y(t,z)D ttZ X£(dz)dt 



< c\\X\ 



for all X 6 D 1 , and c being a constant only depending on Y. If V G Dom(<5), then the 
Skorohod integral of Y, denoted by 

-T [■ 

Y(t,z)N(5z,5t) (18) 



o 

is defined as the unique element in L 2 (P) characterised by 
rT 



E 



X 



o 



Y(t,z)N(6z,5t) 







= E 






.JO Jf 



Y(t,z)D t , z Xl{dz)dt 



(19) 



for all X e P 1 ' 2 . 

The Skorohod integral can be viewed as a creation operator on the chaos: it can be shown 
that Y S Dom(<5) if and only if 



5> + l)!|/ n | 



(20) 



n=0 



where f n G L 2 (([0,T] x M 



is the symmetrization of the nth chaos kernel function 



f n (-,t,z) of Y. The Skorohod integral becomes 

[ T [ Y(t, z)N(Sz, St) = I n +1 (fn) ■ 
Jo Jr n=0 



(21) 



We see that the definitions of the Malliavin derivative and Skorohod integration is completely 
analogous in the Wiener and pure-jump Levy cases. 

We now move on to some properties of the Skorohod integral, which turn out to be slightly 
different from the Wiener case. We start with the "fundamental theorem of calculus" : 

Proposition 3. Suppose that Y £ Dom(5) HID 1 ' 2 . If the random field (s,y) i— > Dt lZ Y(s,y) is 
Skorohod inteyrable, for almost every (t,z), then 



D 



t.z 



(22) 



Y(s,y)N(5y,5s) = Y(t,z)+ / D t>z Y{s, y)N{6y, 8s) . 

10 JR Jo JR 

The next result, which is an integration by parts formula, is particularly useful in our 
stochastic integral definition for VMCV processes. 

Proposition 4 (Integration by parts formula). Suppose that Y E Dom(5) and X £ D 1 ' 2 . If 
the random field Y(t, z)(X + D t , z X) is Skorohod integrable, then 



X 



Y(t,z)N(8z,5t) 



Y{t,z)(X + D t , z X)N(5z,5t) 



+ / / Y(t,z)D t , z X£(dz)dt. (23) 
Jo JRo 

We observe that this relationship is not completely analogous to the Gaussian case, as we 
have the additional term Y(t, z)Dt tZ X inside the Skorohod integral. 
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3 Volatility modulated Levy— driven Volterra processes 



Let us now define the class of volatility modulated Levy-driven Volterra (VA4CV ) processes. 
As before, we denote by L = (L t )t>o a Levy process. A VM.CV process X = (X(t)) t >o, is 
defined by 

X(t)= I g(t,s)a{s)dL(s), t > 0, (24) 
Jo 

where g is a real-valued measurable function defined on the space 

]R2 : = {(t,s) G M. 2 \t > s > 

0} and a = (a(t)) t >o is an ^-adapted, predictable stochastic process. Note that generally 
g(t, s) is not defined for t = s. 

We always assume that the integrand s t— > g(t, s)cr(s) satisfies the integrability conditions 
© a.s. for all t £ M, making X well-defined and 5~adapted. 

Note that in financial modelling, a is usually referred to as the stochastic volatility, while 
in applications to turbulence, it is referred to as the stochastic intermittency. In some appli- 
cations, the stochastic processes a and L are assumed to be independent. 



3.1 Examples of the kernel function 



Let us study some relevant examples of VAiCV processes, which are often used in applications. 

An important type of kernel function is the so-called shift-kernel, which is defined as 
g(t, s) = g(t — s) for all t > s > 0. This choice is motivated from the class of Levy semista- 
tionary (CSS) processes. In the Brownian case such processes have been used for modelling 
turbulence, see lBarndorff-Nielsen &; Schmiegell (12009) . and the ex tension to more general Levy 
processes has been introduced in iBarndorff-Nielsen et al. (l2O10h in the context of modelling 
energy spot prices. 

Example 4. The CSS processes encompass many existing dynamical models. For example, 
by choosing g(t,s) = exp(— a(t — s)) for s < t with a > and a(s) = 1, we recover the 
solution of the Ornstein-Uhlenbeck (OU) process 

dX(t) = -aX(t) dt + dL(t) . 

Example 5. A generalization of OU-processes and still within the class of CSS processes, is 
the class of CARMAfp, pp -processes . A CA RMA process is the continuous-time analogue of 
an ARM A time series, see Brockwell d2001ol fth for details, and it is defined as follows Suppose 
that for p, q € No with p > q, X(t) = b'V(t), where bslP and V(t) is a p-dimensional OU 
process of the form 

dV(t) = AV(t)dt + e p dL(t) , (25) 



with 







-OLr, 



-a 



P-i 



Here we use the notation I p _i for the (p — 1) x (p — l)-identity matrix, e p is the pth coordinate 
vector (where the first p—1 entries are zero and the pth entr y is 1) and h' = \b n, b\ 



'p-U 



is the transpose of b, with b„ = 1 and bj = for q < j < p. iBrockwelll hood ) shows that if 



all the eigenvalues of A have negative real parts, then V(i) defined as 



V(t)= / e A ^ s) e p dL(s) 

J — oo 
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is the (strictly) stationary solution of (|25p . Furthermore, 

X(t) = b'V(t) = f b'e^-^epdLis) , 

J—oo 



(26) 



is a CARMA(p, q) process. Hence, a CARMA process is essentially an CSS process with 
constant volatility and kernel function specified as g(t,s) = b'exp(A(£ — s))e p . 

Example 6. In applications to turbulence, one often works with the following function: 

g(t, s) = g{t -s) = {t- s) u ~ l exp(-A(t - s)) , (27) 

with v > | and A > 0, see iBarndorff-Nielsen fc Schmiege 3 (|2008h . This choice leads to a 
well-defined CSS process, which is stationary (provided the stochastic volatility component 
is stationary). However, notice that g(-) has a singularity at zero when v € 1). 

Example 7. Important examples within our modelling framework are also fractional Brow- 
nian motion and fractional Levy processes. To see this, suppose that a = 1. Choosing 



g(t,s)=c(H)(t-s) H - 1 / 2 + c(H) 
with 

c{H)- 



1 



H 



\ u _ s f-3/2 ^ _ (s / n )V2-^ du ) (28) 



2flT(| - if) 



/ r(ij + i)r(2 - 2H) ' 

and i/ G (0, 1), we recover fractional Brownian motion (see Alos et al. ( 200ll ). page 798): 



X(t) = [ g(t,s)dB(s), t > 0, 
Jo 



where B is a standard Brownian m otion. Recall th e following alternative representation of 
fractional Brownian motion (cf. e.g. Nualart ( 2006 J )): 



where 



and 



g(t,s) 



Ci(H) 



X(t) 
(t- s 



g{t,s)dB{s 
H-l/2 , \H-l/2 



x + := max{x, 0}, 

1/2 



(29) 
(30) 



Replacing B in (f29]l by a zero mean square-integrable two-sided pure-jump Levy process L(t) 
and a ssuming H 6 (1/2, 1) we get a so-called fractional Levy process as defined bv lMarquardt 
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3.2 Examples of the stochastic volatility /intermittency process 

So far, we only mentioned that we allow for stochastic volatility /intermittency. Let us briefly 
point out relevant specifications of such processes. We often model a 2 as a VM.CV process 
itself, i.e. let 

t 

i(t,s)dL a (s), (31) 



o 

where i is a real-valued measurable function defined on IR> which is integrable with respect 
to the Levy subordinator L a . Moreover, if one wants to ensure that the volatility process is 
stationary, one can work with the specification 

oi= [ i(t - s)dL a (s), (32) 



where L a is a two-sided Levy process and the integration is defined as in iRaiput fc Rosiiiski 



(1989). A simple example which involves dependence between a and L is the following one. 



Example 8. Let L denote a two-sided Levy process with inverse Gaussian (IG) marginal 
distribution. Further choose L a = L and i(t, s) = exp(— (3(t — s)) for /3 > 0. Then a 2 is 
an Ornstein-Uhlenbeck process driven by an inverse Gaussian Levy process, i.e. an IG-OU 
process. 

3.3 Semimartingale conditions 

It is important to note that generally VMCV processes are not semimartingales. E.g. consider 
the case of Example [6l which is relevant in the context of modelling turbulence; there the 
kernel function was chosen as g(t, s) = (t — s) 1 " -1 exp(— X(t — s)), with v > \ and A > 0. When 
v € (i, 1) or v G (1, |], the corresponding VM.CV process is not a semimartingale, and these 



cases are of primary interest in the turbulence context, see iBarndorff- Nielsen &; Schmiegel 
(120091 ). 



However, as long as g is non-singular and satisfies a smoothness condition, we can show 
that X is a semima rtingale. Indee d, we have the follow ing result, w hich is closely related to 
e.g. Corollary 4.9 in IProtterl (jl985l l and Theorem 4.6 in basse! (120091 ): 



Proposition 5. Suppose X is a VAiCV process defined as in (|24j) . and assume that t \— > a(t) 
is locally bounded pathwise a.s. Suppose g(t, s) is defined for all < s < t and that there is a 
bi-measurable function tfi(t, s) such that 

g(t, s) = g(s, s) + / 4>(v,s)dv forall0<s<t, (33) 

J s 

where f^igis, s)) 2 < oo for all t > and 

"t fU 

4>(u, s)) 2 ds du < oo, for all t > 0. (34) 



'o Jo 

Then, (X(t))t>o is a semimartingale with decomposition 



r-t 

dX(t) = g{t,t)a(t)dL{t) + / <p(t, s)a(s) dL(s) dt , (35) 

o 
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that is 



X 



(*)= f 9{ 
J o 



s, s)cr(s) dL(s) + 



j o 



^(v, s)a(s) dL(s) dv, 



t > 0. 



(36) 



Note that (|33p essentially is equivalent to the assumption that g(s, s) exists and t i-> s) 
is differentiable for t > s with <f>(t, s) = |f (i, s). 



Proof. We have 



X(t) 



g(t,s)a(s)dL(s) 



g(s, s)a(s) dL(s) + 



g(s, s)a(s) dL(s) + 
4>(v, s)a(s) dL(s 



<fi(v, s) dv 



a(s) dL(s) 



dv, 



where we have applied Fubini's theorem, see lProtterl (120051 . Theorem IV. 65), which is justified 
by the assumptions on <fi and a. □ 

See also Basse fc Pedersen ( 20091 ) and Basse-O'Connor ( 2010l ) for related results. 

Example 9. Consider the example from turbulence where we choose g(t, s) = (t—s) u ~ 1 e~ x ( t ~ s * > 
with A > and v > ^. Note that assuming that v > | implies that g(s,s) = 0. The in- 
tegrability conditions ([33]) and (f3i|) are then satisfied with 4>(t,s) = ^j-(t,s) for t > s. Also 
note that t i— > g(t, s) is differentiable at t = s when v > 2. According to Proposition [U the 
semimartingal e cases are v = 1 and v > %. In f act, when v > ^ the process is even of bounded 
variation. See Barndorff-Nielsen fc Schmiegel ( 20091 ). Barndorff-Nielsen et al. ( 2011 ) for more 
details. 



4 Stochastic integration 

In this section we are concerned with stochastic integration with respect to a VM.CV process 
X, that is, to define 

f Y(s)dX(s), t>0, (37) 
Jo 

for a class of stochastic processes Y(s) , s E R. The integrands Y may be VMCV processes 
themselves. We are concerned with having a definition which may facilitate the integration 
with respect to possibly singular VA4CV processes, in the sense of kernel functions g(t, s) 
which may not be defined for t = s. 

Since at present we only have a theory of integration with respect to Brownian and pure- 
jump Levy processes separately we consider two cases: First, we study the case when X is a 
volatility modulated Brownian Volterra (VA4BV) process, 

X(t)= [ g(t,s)a(s)dB(s), (38) 
Jo 

where B is a standard Brownian motion. Next, we treat the pure-jump VM.CV case 

X(t) = f g(t,s)a(s)dL^(s), (39) 
Jo 

where is a zero-mean square integrable pure-jump Levy process. 
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4.1 The VMBV case 

We start with the case of a volatility modulated Brownian driven Volterra process. Through- 
out this Section, we will be working under the following assumption. 

Assumption A The process B = (B(t))t>o is a standard Brownian motion. For t > let Tt 
be the cr-field generated by B(s) with s <t. The process (a(t))t>o is (J r t) t >o-predictable 
and ^ 

E( f (g(t,s)a(s)) 2 ds) < oo, (40) 



for all t > 0. Finally, J- is the cr-field generated by B. 

When Assumption A is satisfied X(t) defined as in (|38p is a well-defined square-integrable 
zero mean VMBV process. 

In order to motivate our definition of stochastic integration, we start out with a heuristic 
derivation. First, supposing t i— > X(t) and s i— > Y{s) are differentiable, integration by parts 
yields 



! Y{s)dX{s) = fv{ 
Jo Jo 



dX{u) t /•* dY(u) 



du = {Y(s)X(s)Y s=0 - / —^X{u) du 



du Jo du 

Y{t)X{t)- f'^M. [ U g( u ,s)a(s)dB{s)du. 







du 

In the integral on the right-hand side, we use Proposition [2] to get 



f ' Y{s)dX(s) = Y(t)X(t)- [ f U ^^-g(u,s)a(s)SB(s)du 
Jo Jo Jo du 

1 f u (dY(u)} 
J Ds \ du j d(u, s)cr(s) ds du . 

Here, 5B denotes stochastic integration with respect to Brownian motion in the Skorohod 
sense, while D s is the Malliavin derivative. Applying a stochastic Fubini theorem, it holds 
that 

" Y{s) dX(s) = Y(t)X(t) - I ( [\(u,s)^^du S \ a(s)6B(s) 



o 



f ( f dY(u) ) 

J D s<^J g{u, s) du du j a(s) ds . 



Here, we used that the Malliavin derivative is linear. On the first term on the right-hand 
side we again apply Proposition [2] to reach 

Y(t)X(t) = Y(t) f g(t,s)a(s)dB(s) 
Jo 

Y(t)g(t,s)a(s)5B(s)+ [ D s {Y(t)} g(t, s)a(s) ds . 



Thus, collecting terms we reach the expression 
f-t ft ( ft 



j Y(s)dX(s) = J\Y(t)g(t,s) - J g( u , s )^^du ] ja(s)SB(s 
+ J D s tY(t)g(t,8)- j\(u, 



— -j— du \ a{s) ds . 
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Observe that after doing an integration by parts, we find 
Y{t)g(t,s)- J g{u,s)^^du 

dg(u,s) 



Y(s)g(s,s)+ [ Y(u 

J s 



du 



du 



= Y(s)g(t, s) + J (Y(u) - Y(s)) du . 

Here we have assumed that g is differentiable with respect to its first argument. Note that in 
the last equality we have reorganized the integral such that the term g(s, s) is not appearing, 
thus ensuring that we can include integration for VAiCV processes which may be singular at 
g(s,s). 

The above derivation motivates the introduction of the operator 

K g (h)(t,s) = h(s)g(t,s)+ [ (h(u)-h(s))g(du,s), (41) 

J s 

whenever integration with respect to g(du, s) makes sense as a Lebesgue-Stieltjes integral. 
That is, for any s > 0, the mapping t i— > g(t, s) is of bounded variation for t in any b ounded in 



terval [u, v] with u > s. We remark that precisely the same operator K, g appears in lAlos et al 



(|200ll . p. 770) in their definition of stochastic integration for Gaussian processes. In the 
par ticular case of fracti onal Brownian motion this operator has also been very useful, cf. 
e.g. Biagini et al. ( 20081 ). Note that the definition does not require the existence of g(s,s). 



Moreover, the mapping 1 1— > g(t, s) can be of unbounded variation on (s, v ) for any v > s. 

Example 10. Note that the example of key relevance in modelling turbulence is the choice 
of the kernel function given in equation (|27p . The above operator can handle this case for any 
v > 1/2. To see that note that t i-4 g(t,s) is of bounded variation on [u, v] but unbounded 
variation for t G (s, v ) for s < u < v if v G (|, 1); further g(s, s) = 1 and t h-> g(t, s) is of finite 
variation on any interval if v = 1, and g(s, s) = and t h- ¥ g(t, s) is of bounded variation on 
any interval if v > 1. In the case of a fractional Brownian motion where g is given by (I28p or 
pOj) . g(s,s) does not exist and t h-> g(t,s) is of unbounded variation on (s,v) for any v > 
when H € (0, ^), g(s, s) = 1 and t i— > g(,s, i) is of bounded variation on compacts when H = \ 
and g(s, s) = and i i-> g(,s, t) is of bounded variation on compacts when H > ^. 

If g(s,s) < oo, we can redefine JC g (h) as long asm-> /i(n) is integrable with respect to 
g(du, s) in the Lebesgue-Stieltjes sense as 

JC g (h){t,s) = h(s)g(s,s)+ [ h{u)g(du,s). (42) 

is 

There is also the case of g(u, s) being absolutely continuous with respect to Lebesgue measure, 
yielding 

dg 



K g (h)(t, s) = h(s)g(t, s) + J (h(u) - Hs))-^(u, s) du (43) 

with dg/du being the Radon-Nikodym derivative of g with respect to the first variable. 
Based on the above heuristic arguments we are now ready to define J" * Y(s) dX(s). 
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Definition 2. Let Assumption A be satisfied and let X be given by (|38p . Assume that for 
s € M+ the mapping u i— >■ s) is of bounded variation on [u, t>] for &110<s<u<v< oo. 
Fix t > 0. 

We say that a process 

s h. y( s ) , 

for s S [0, i] belongs to l x (0,t) if the following conditions are satisfied: 

1. For s € [0,t] the process (Y(u) — Y(s)) ue ^ s ^ is integrable with respect to g(du,s) a.s. 

2. The mapping 

s^JC g (Y)(t,s)a(s)I m (s), 
is Skorohod integrable with respect to B. 

3. JC g (Y)(t, s) is Malliavin differentiable with respect to D s for s E [0,i], with 

s^D s {K g (Y)(t,s)}a(s) 
being Lebesgue integrable on [0,t]. 

Suppose that s h-> ^(s) for s 6 [0, i] belongs to X x (0,t). Then the stochastic integral of Y 
with respect to X is defined as 

I Y{s) dX{s) = I lC g (Y)(t, s)a{s) 5B(s) + / D s {K g {Y){t, s)}a(s) ds . (44) 
Jo Jo Jo 

The integral on the right-hand side of (1441) is defined in the sense of Subsection 12.11 where 
as always we use the convention J* := J* R I[o,tl- 

There is no conditions of adaptedness on the stochastic process Y(s). Our definition 
gives an anticipative stochastic integral. Note that the definition does not make use of g(t,t) 
explicitly, and thus we can include kernel functions which are singular, such as the choice 
of kernel function often used in modelling turbulence, see (|27|) . as kernels in the Brownian 
semistationary process. 

Note h e re tha t the integral we define in (|44p generally differs from the integral studied by 
Alos et al.1 (|200lh since we have the additional integral with respect to the Lebesgue measure. 



4.2 The pure— jump VAiCV case 

Now we turn our attention to pure-jump VMCV processes. In what comes, we will work 
under the following assumption. 

Assumption B Let T < oo denote a finite time horizon. The process = (£^(t))tG[o,Tl 
is a square integrable pure-jump Levy process with zero mean and, for < t < T, Ft 
is the (j-field generated by L^ p \s) with s < t. The process (<7(i))te[o,T] is {^Ft)t&[o,T\- 
predictable and (l4T)j) is satisfied for all t € [0, T]. Finally, T = Ft- 

Let us work under Assumption B and let (A"(t)) tg r 0i Ti be the pure-jump VMCV process 
defined by (j39|) . By decomposing as in (|15j) . X is given by 

X(t)= [ g(t,s)a(s)dL ( P\s) = [ [ zg(t, s)a(s) N(dz,ds). 
Jo Jo it 
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Our goal is the same as in the VA4BV case, to define 

t 

Y(s)dX(s), 



o 

for a class of integrands Y(s). The extension to the Le vy case rests on the Malliavin Calculus 
for Levy processes. We refer to Di Nunno et al. ( 20051 ) for the development of such a calculus, 
cf. also Subsection 12.21 The main ingredients of interest to us are the extension of Skorohod 
integration, the Malliavin derivative and the integration by parts formula for constants. 
We start first with a heuristic derivation following the VM.BV case: We find 

Y(s)dX(s)=Y(t) /"* g(t,s)a(s)dL ( Phs) - [* lUM. g( u , s )a(s) dL& (s) du . 
Jo Jo du J 

By applying the integration by parts formula in Proposition HI we find for the first term 

Y(t) [ [ zg(t,s)a(s) N(dz,ds) = [ [ zg(t, s)a(s) N(dz,ds) 

o Jr 

+ / f zg(t, s) (Y(t) + D s , z {Y(t)}) a{s) N(5z, ds) 
zg{t,s)D S:Z {Y(t)}a{s)£{dz)ds. 



o 



Here, D S)Z is the Malliavin derivative with respect to the Poisson measure N, and N(5s, dz) 
denotes the Skorohod integral with respect to the compensated Poisson random measure. 
Note the additional term D StZ {Y{t)} in the Skorohod integral which is not present in the 
corresponding term in the VMBV case. Similarly, we find for the second term 



dY(u) 
du 



zg(u, s)a(s) N(dz, ds) du 
o Jo Jr 



du ' i du 

+ Jo So L z9{u ' s)Ds,z { } a{s)e{dz) ds du 

r \ g (u,s)(^^ + D s ,J^-X\ dua(s)N(5z,5s) 

+ Si S SI ' zg ^ s)Ds ' z \^r\ du£ ( dz ) a ( s ) ds ' 

where we have used the stochastic Fubini theorem in the last equality. Next, by the linearity 
of the Malliavin derivative, we find 



dY{u) 
du 



zg(u, s)a(s) N(dz, ds) du 

( f* , ,dY(u) , _ f f* , ,dY(u) , 
z[ / g(u,s) , v ; du + D S)Z { I g(u,s) , v ' du\ a(s) N(5z,5s) 



du ' | ./« du 

•t 



+ f R zD s,Aj s 9 ^ s ^^~ du ] t( dz M s ) ds > 



17 



As a last step, by again doing an integration by parts as in the VMBV case, and collecting 
terms, we reach 

Y{s)dX(s)= f* f z(JC g {Y){t,s) + D Sj2 {]C g (Y)(t,s)})a{s)N{Sz,ds) 

Jo Jr t (45) 

+ / / z(D SiZ {lC g (Y)(t, S )}) e(dz)a(s)ds. 
Jo Jr 

This relation involves a Malliavin derivative for Levy processes and a Skorohod integral for 
Poisson random measures, along with the operator K, g . 

Next, we introduce the class of integrands T x (0,t) for stochastic integration with respect 
to pure-jump VJvlCV processes. 



Definition 3. Let Assumption B be satisfied and let X be given by f|39jl . Assume that for 
s E [O.T] the mapping g(u,s), defined for u G (s,T\, is of bounded variation. Fix t g]0,T]. 
A process 

s H. Y(s) , 

for s G [0,t] belongs to I (0, t) if the following conditions are satisfied: 

1. For s G [0,t] the process (Y(u) — Y(s)) u( z^ s ^ is integrable with respect to g(du,s) a.s. 

2. JCg(Y)(t,s) is Malliavin differentiable with respect to D s>2 for (s,z) G [0,t] x R, with 

(s,z)^zD s , z {K g (Y)(t,s)}a(s), 

being t{dz) cis-integrable. 

3. The mapping 

(s, z) ' y z (K g (Y)(t, s) + D s>z {K g (Y)(t, s)}) a(s)I m (s), 

is Skorohod integrable on [0,T] x R with respect to N(ds,dz). 

Suppose that s h-> Y(s) for s G [0,t] belongs to l x (0,t). Then the stochastic integral of Y 
with respect to X is defined as 

Y(s)dX(s)= [ [ z(JC g (Y)(t,s) + D s , z {JC g (Y)(t,s)})a(s)N(5z,5s) 
o Jo Jr 

+ f [ zD s , z {lC g (Y)(t,s)} l{dz)a{s)ds . 
Jo Jr 

The first integral on the right-hand side is defined in the sense of Subsection 12.21 where 
Jo := Jo f [0,t]- 

4.3 Further remarks 

In order to verify conditions (1) in Definition [2] and [3] we must study existence of K, g {h). 
Recall that ^ 

Kg(h)(t, s) = h(s)g(t, s) + / (h(u) - h(s))g(du, s). 
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Thus, K,g(h)(t,s) exists if e.g. h(t) = zl{ a <t<&} for some constants z G R and a < b. If 
the measure g(du,s) is finite on any finite interval then fC g (h)(t,s) exists if h is bounded 
on compacts and in particular if u i — h(u) is cadlag or caglad. This covers the turbulence 
example in equation (|27p with v > 1 and fractional Brownian motion with H > 1/2. In 
general, however, some kind of smoothness of h(t) near t = s is needed to ensure existence 
of fC g (h)(t,s). For example, assume u i-> g(u,s) is C 1 with ^g{u,s) ~ (u — s) a as m | s 
where a > —2. If there is a constant X such that \h(u) — h(s)\ < K(u — s) for all u > s, then 
K.g(h)(t,s) exists. This covers the turbulence example (equation (|27|) ). where a = v — 2, as 
well as fractional Brownian motion, equations (|28p or (|30p . where a = if — 3/2. 
The following Lemma gives useful information in the CSS case. 

Lemma 1. Suppose g is a shift-kernel, i.e. g(t,s) := g{t — s). If f is Lebesgue-Stieltjes 
integrable on [s,t] with respect to g(du,s), then f(s + •) is Lebesgue-Stieltjes integrable on 
[0, t — s] with respect to g(du) and 

1 f(u) g(du, s) = I 3 f(s + u) g(du) . (46) 
J o 

Proof. Note that if E G B([s,t}), then F := E - s G B([0,t - s)), where B(A) is the Borel 
cr-algebra on A, an interval of R. Since for any (a, 6] G £>([s, t]) we have that (a — s, b — s] G 
B([0,t- s\), it holds that 

c/((a, 6], s) = 5(6 - s) - g(a - s) = g((a - s, b - s]) . 

implying translation invariance of the Lebesgue-Stieltjes measures g(d-,s) and g(d-). Hence, 
by definition of the Lebesgue integral the Lemma holds. □ 

The Lemma gives us the form of fC g in the case of g being a shift-kernel: 

rt— s 

K g (h)(t,s) = h(s)g(t-s)+ (h(s + u) - h(s)) g(du) . (47) 

Jo 

Note that JC g (h) does not become a shift kernel. 



5 Some properties and applications 

In this section we will either work under Assumption A, in which case X is given by (|38p . 
or under Assumption B, in which case X is given by (|39p . Let L = B when working under 
Assumption A and L = when working under Assumption B. 

5.1 Fundamental properties 

Let < t\ < t2 and assume s h-> Y(s) is in I x (0, t\) (resp. ml x (0, t\)). Then s i— > Y(s)I{ s<tl } 
is in T x (0, £2) (resp. in I x (0,t2)), and we have the usual rule 

[ 2 Y(s)I {s < tl} dX(s) = I' Y{s) dX{s). 
Jo Jo 
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The definition of a stochastic integral is linear: If Y and Z are two processes in l x (0,t) 
(resp. in X (0, t)) and a, b two constants, then we easily see by linearity of the Malliavin 
derivatives, Skorohod integration and Lebesgue integration that 

/ (aY(s) + bZ(s))dX(s) = a [ Y(s)dX(s) + b f Z{s)dX(s). 
Jo Jo Jo 

This in particular implies that if < u < v < t and s i— y Y(s) is in I x (0, u) and in I x (0, v) 
(resp. in l x (0,u) and in l x (0,v)), then s i->- Y(s)I^ u<s < v y is in X x (0, t) (resp. in T x (0, i)) 
and ^ 

Y(s)I {u<s < v} dX(s)= [ l Y(s)dX(s)- [ U Y( S )dX(s). 



Since it is obvious that IC g (l)(t, s) = g(t, s) for s < t and D s K g (Y)(t, s) = U SjZ K, g (Y)(t, s) = 
when Y" is deterministic, it follows that 

/ ldX(s) = X(t) and, more generally, / ^-{u<s<v} 

dX(s) = X(v) - X(u). (48) 

Jo Jo 
The following integration-by-parts formulae hold: 

Proposition 6. (1) Let Assumption A be satisfied and X be the VMBV process given by 
(|38D . Let s i-> Y(s) be an element in l x (0,t). Suppose Z is a bounded random variable such 
that the process s i— > ZY(s) G T (0, t). T/ien i£ ZioWs i/iat 



/ ZY(s)dX(s) = Z [ Y(s)dX(s). 
Jo Jo 



(49) 



^ Lei Assumption B be satisfied and X be the VMCV process given by (j39j) . Let s \— > Y(s) 
be an element in i x (0,t). Suppose Z is a bounded random variable such that the process 
s ^ ZY(s) € l x (0,t). Then it holds that 



f ZY(s)dX(s) = Z f Y(s)dX(s) 
Jo Jo 



[ [ zD s . z {K g (Y)(t, a)} D SjZ {Z} a{s) (n(5z, Ss) + £{dz) ds 
Jo Jr v 



(50) 



Proof. Without loss of generality, we suppose for notational simplicity that a(s) = 1. Note 
that 

K g (ZY)(t,8) = ZlC g (Y)(t,8). (51) 
Proof of (1). By the product rule for the Malliavin derivative 

D s {ZK g (Y)(t, s)} = D s {Z} K g (Y)(t, s) + ZD S {K g (Y)(t, s)} . (52) 

Hence, by (f5"Tj) we have 

f ZY{s)dX(s)= [ ZIC g (Y)(t,s)6B(s) 
Jo Jo 

+ / ZD s {K g {Y){t,s)} + D s {Z}K g {Y){t,s)ds. 
Jo 
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Applying Proposition [2] this gives 

f ZY(a)dX(a) = z( f IC g {Y){t,s)5B(s)+ [ D s {lC g (Y)(t, S )} ds 
Jo \Jo Jo 

= Z I Y(s)dX(s), 
Jo 

which is (1). 

Proof of (2). Bv lDi Nunno et al.1 (|2005l . Lemma 3.11) we have (|52p with D s replaced with 
D s z . Hence by this equation and (j5"Tj) , 



f ZY{s)dX(s) = I I zK g (Y)(t,s)(Z + D s>z {Z}) N(5z,6a) 
Jo Jo Jr 

+ 11 zD StZ {JC g (Y)(t, a)} (Z + D s>z {Z}) N(5z, 6s) 
Jo Jr 

-II zD s>z {K g (Y)(t, a)} D s . z {Z} N(Sz, 5s) 
Jo Jr 

+ f f z (ZD S:Z {IC g (Y)(t, s)} + D S:Z {Z} K g {Y){t, a)) £(dz) ds 
Jo Jr 



Invoking Proposition H] gives 

/ ZY(s)dX(s) = Z [ [ zKg{Y)(t,s)N{Sz,6a) - [ [ zfC g {Y){t, s)D 3>z {Z} £(dz) ds 
Jo Jo Jr Jo Jr 

+ Z [ [ zD SyZ {K g (Y){t,s)} N(5z,5s) 
Jo Jr 

zD SyZ {K g (Y)(t, a)} D s z {Z} £(dz) ds 

-II zD s , z {JC g (Y){t,s)}D S:Z {Z} N(8z,5s) 
Jo Jr 

+ I I z (ZD SyZ {K g (Y)(t, a)} + D SyZ {Z} JC g (Y)(t, a)) £(dz) ds . 
Jo Jr 

Hence, the result follows by collecting terms. □ 

We note that (|49p is a classical rule of calculation, which, by Proposition[6](2), is generally 
only satisfied in the Gaussian case. But (|49p is also satisfied in the pure-jump case when the 
integrand Y is deterministic since we then have D s ^ z {JC g (Y)(t, a)} = 0. 

Let us consider the integration of "simple" integrands. Suppose that < io < h < £2 < 
• • • < t n < t is a partition of the interval [0, t] and let Zi be bounded random variables such 
that s 1 y Zilr t . <s<t 1 belongs to l x (0,t) respectively to l x (0,t). Then we call 

n-l 

y ( s ) = Ys ZiI {u<s<k +1 } 

i=0 

a simple integrand. Since s 1— > \ti <s<t i+1 } is deterministic it follows from Proposition El 
equation (|48p and linearity of the integral that 

n-l 

/ Y(s)dX(s) = Y j Z i AX(t i ), (53) 
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with AX(ti) = X(ti + i) — X{ti). Thus we see that the integral coincides with the simple 
integral on simple (non- adapted) integrands. We remark that thi s property gener ally only 
holds for deterministic simple integrands for the integral defined in (jAlos et al.ll200ll . p. 773). 

We next show that our integral shares the localization property, which is also one of the 
basic properties of the ltd integral. 

Proposition 7. Suppose that Y{s) = 0, a.e. s < t, a.s. Then, under Assumptions A or B, 

[ Y{s)dX(s) =0,a.s., 
Jo 

Proof. Since Y(s) = for a.e. s <t, a.s., we have that 

(Y(u)-Y(s))g(du,s) =0, 

for a.e. s < t a.s. Thus, K, g (Y)(t, s)a(s) = for a.e. s < t, a.s. Then, by Nualart ( 2006^ . 
Proposition 1.3.6 ) we find that 

lC g (Y)(t,s)a(s)5B(s) = 0, a.s. 

Furthermore, we have D s JC g (Y)(t, s) = and D SjZ JC g (Y)(t, s) = 0, except on sets of measure 
zero. Hence, using the property of Lebesgue integrations, we find that 

J (D s K g (Y)(t,s)+ J^zD s>z K g (Y)(t, S )e{dz)a(s)\ ds = 0, a.s. 

Finally, by making use of Corollary 3.9 in lDi Nunno et al.l (|2005l ). we find that 

*t 



□ 



o 



z {fC g (Y)(t, s) + D S)Z {JC g (Y)(t, s)}) a(s) N(5z, 8s) = 0. 
Hence, the result follows. 



5.2 Applications 

5.2.1 Deterministic integrands 

In this subsection we focus on deterministic integrands on the form Y(s) = h(t,s), with 
s i y h(t, s) being a measurable function on [0, t] such that JC g (h)(t, s) is well-defined (that is, 
h(t,u) — h(t,s) is integrable with respect to g(du,s)). Since we know that in this case the 
Malliavin derivatives in (j44j) and (|45p become zero, we have that 

t rt 
h(t,s)dX(s)= / K g {h){t,s)a{s)dL{s), 
Jo 

where we recall that L = B or L = depending on which Assumption we are working 
under. Thus, the integral of a deterministic function with respect to a VM.CV process gives 
again a VA4CV process. 
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If we consider shift-kernel functions h,g, i.e., h(t, s) = hit — s),g(t, s) = git — s), we find 
that JC g (h)(t, s) depends on s and t only through (t — s), since from (|46D 

rt—s 

K g (h)(t, s) = h(t - s)g(t - s) + / (h(t - s - u) - h(t - s)) g{du) . 

Jo 

Let us look at a simple example: Consider h(u) = exp(— au), g{u) = exp(—/3u) for a, j3 > 
with a 7^ j3. Then, after a straightforward integration, 



and we find 

f t e-^- s Ux(s)= f\ 

o Jo {a- 13 a-P 



1 e - a{t ~ s) dX{s) = f { -^- e - a ^ ^L e -/3(*-s) \ a (s) dL{s) 



If a = P, we trace back and find 



e -a(t~s) dx ^ = / e -«(^)(l_ a (t_ s )) (J ( s )dL(s). 
JO 

Remark. This type of integral kernel appears in connection with stochastic volatility mod- 
elling for commodity prices based on the Schwartz dynamics and the Barndorff-Nielsen and 
Shephard model, see iBenthl (|201ll ) . 

Moreover, as the next Proposition shows, the integral of h with respect to X is also the 
solution of an V.M£V-driven Ornstein-Uhlenbeck process: 

Proposition 8. Assume that for all s € M + the mapping (s, oo) 3 u \— > e~ a ( u - s ) - 1 b 
integrable with respect to g(du,s), and 

lim(l - exp(-a(i - s)))g(t, s) = . (54) 
Assume in addition that t h-> a(t) is pathwise locally bounded a.s. Finally, assume that 



t rt 



o Ju 



(g(s, u)) 2 ds du < oo . (55) 



Then, s t— > e~ a ^~ s ^ is in I x (0,t) for all t > when Assumption A is satisfied (or in X x (0,t) 
for all t <T when Assumption B is satisfied) and Y(t) = J * exp(— a(t — s)) dX(s) solves the 
VA4CV -driven Ornstein-Uhlenbeck process 

Y{t) = -a[ Y(s)ds + X(t). (56) 
Jo 

Proof. By Assumption, (1) in Definition [2] (resp. in Definition [3]) is satisfied and since the 
integrand is deterministic this implies that s *— > e -Q ( i-s ) is in l x (0,t) (resp. in l x (0,t) for 
all t < T). 
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Denote for simplicity h(u) = exp(-au), and we deduce from the definition of integration 
with respect to X that 

ft ft ft fS 

Y(t) + a Y(s)ds= / h{t-u)dX{u) + a / h(s - u) dX(u) ds 
Jo Jo Jo Jo 

t ft fS 

JC g (h)(t,u)a(u) dL(u) + a / / K. g (h)(s,u)a(u) dL(u) ds, 
o Jo Jo 

where we recall that L is either B or . Appealing to the stochastic Fubini Theorem cf. 
details below we find 



Y{t) + a J Y(s)ds = J (jC 



g (h)(t,u)+a I K. g (h)(s,u) ds \ a(u) dL(u). 



t r t 



Now, we have that 

K g Qi)(t, s ) + a f s Kg( h )(u, s) du = e~ a{t - s) g(t, s) + (^ e ~ a{t ~ u) - e~ a{t ~ s) ^j g{du, s) 

+ aj* (e- a{u - s) g(u,s) + £ (e~ a ^ - e - Q ("- s )) g(dv,s) \ du . 
But, by the Fubini theorem, we find 

01 f. C { e ~ a(U ~ V) ~ e" a( "" s) ) g(dv, s) du 

^-*(u-v) _ e - a (u-s)\ dug ( dViS ) 
^ e -a(t-v) _ e -a(t- S )^ gfdv ^ s) + J ^ _ e -a{v-s)^ g ^ g j _ 

After doing an integration-by-parts and applying the limit condition on g, we get 

(l _ e -«(«--)^ g (dv, s) = (1 - e- a ^)g{t, s) - a J e^ 1 *-^ g(v, s) dv . (57) 
This gives that 

K g (h)(t, s) + a J ICg(h)(u, s) du = g(t, s), 

and the proposition is proved except for the verification of the Fubini theorem, which is 
relegated to the following Lemma. □ 

T o justify the u se of stochastic Fubini above we have to verify the following condition 
from lProtterl (|200.4 Theorem IV.65). 



Lemma 2. Using the same notation as in the proof of Proposition^ 



t r t 



o Ju 



(ICg(h)(s,u)) 2 ds a 2 (u) du < oo, a.s. (58) 
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Proof. Since u i— > o~(u) is locally bounded we may and do assume that o~(u) = 1. By definition 
of K g (h) and (|57|l. 



K. g (h)(s, u) = h(u)g(s, u) + / (h(v) — h(u)) g(dv,u) 

J u 

= h(u)g(s, u) + e- QU / (e"^-^ - 1) 5 (cfo, «) 

J u 

= h{u)g(s,u) - h(u){(l - e~ a( - s ~ u ^)g(s,u) - a \ h(v - u)g(v,u) dv} 
= h(s)g(s,u) + a / h(v)g(v,u) dv . 



Thus, it suffices to verify (|58p with fC g {h) replaced by the two terms on the right-hand side. 
Since h is bounded we have by ([55]) that 



t i-t 

I (h(s) g (s, u)) 2 ds du < oo . 

'0 Ju 

Since h(v) is bounded by 1 for u>0we have by Jensen's inequality that 

t rt / p8 \ 2 rt rt rs 

I I / h(v)g(v, u) dv J dsdu< / / (g(v,u)) 2 (s — u) dv ds du 

Ju \J u J J J u J u 

rt rt rt rt rt 

<t / / / (g(v,u)) 2 dv ds du < t 2 l j (g(v, u)) 2 dvdu < oo 

J J u J u Jo Ju 

which gives (f58|) . □ 

We notice that a sufficient condition for the limit condition in the proposition to hold is 
that the limit of g(t, s) exists as t X s. 

Example 11. Recall the example from turbulence: 

g(t, s) = g(t -s) = (t- s)"- 1 exp(-A(t - s)) 

for v > |. This gives a singularity at s = t when v € 1). However, appealing to the 
L'Hopital rule, it is easily verified that the condition (|54p is satisfied. Hence, the condition on 
the limit behaviour of g is satisfied for all v > \. Similarly, for fractional Brownian motion 
with kernel (|28p or (|30p there is a singularity at t = s when H £ (0, ^) but the same kind of 
argument shows that the limit behaviour of g is satisfied for all H € (0, 1). 

Note that in general, the solution Y(t) of the V.M£V-driven OU process in (|56p is not 
in a shift-kernel-type form. Since g may not be a shift-kernel, we also get that Y{t) is a 
VAiCV process with a kernel function which is not of shift-kernel type. If g is a shift-kernel, 
then also Y(t) will be a (semi) stationary VA4CV (or more precisely an jCSS) process. In 
Brodv Sz Zervosl ( 2002 ) one has suggested Ornstein-Uhlenbeck processes driven by fractional 



Brownian motion as a model for t he dynamics of temperature. Our ana l ysis h ere links their 
results to CSS processes. See also Barndorff-Nielsen Sz Basse-O'Connor ( 201ll ) for a general 
treatment of so-called quasi Ornstein-Uhlenbeck processes and references to the literature in 
this field. 
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5.2.2 A Volterra process as integrand 

It is natural in the present context to consider the question of the characteristics of integrals 
of Volterra processes with respect to Volterra processes. We shall not aim here to give a 
comprehensive discussion of this point, but will restrict to an illustration in the VA4BV case. 
This will be connected to chaos expansions with respect to X and, as we shall see, this "double 
integral" with respect to X gives an element in second chaos, as expected. 

Throughout this section we work under Assumption A and X is given by (|38p . Further 
assume that a it) = 1, i.e., there is no stochastic volatility, and that s *— > g(t, s) is a measurable 
function on [0,t] where g(s,s) is well-defined for all s € [0,t]. Consider the integrand 



Y(s) = / h(s,v)dX(v) , 
Jo 

for some sufficiently regular function h G I x {0, s). We want to compute the integral 

Y(s)dX(s) 



o 

by using our theory for stochastic integration. 
By Definition [2] we find that 

f-t ft 



[ Y(s)dX(s)= [ K g (Y)(t,s)6B(s)+ [ D s {K g (Y)(t, s)} ds , 
Jo Jo Jo 



where 

Y(s) = f S fC g (h)(s,v)dB(v). 
Jo 

First, let us compute JC g (Y)(t, s): By definition we have 
K g {Y)(t t s) = Y(s)g(t, s) + f (Y(u) - Y(s)) g(du, s) 

J s 

K g (h)(s,v)g(t,s)dB(v) 

+ J* ^ ]C g (h)(u, v) dB(v) - J* JC g (h)(s, v) dB(v) \ g(du, s) 

lC g (h)(s,v)g(t,s)dB(v) + / K g (h)(u,v) - K. g (h)(s,v) dB(v) g(du,s) 

Js Jo 

+ / K g (h)(u,v) dB(v) g(du,s) 

J s J s 

K g (h)(s,v)g(t, s) dB(v) + / / K g (h)(u,v) — K. g (h)(s,v) g(du,s) dB(v) 

Jo Js 

+ / / K. g (h)(u,v) g(du, s) dB{v) . 



In the last equality we applied the stochastic Fubini theorem, cf. e.g. iProtterl (|2005l . Theorem 
IV. 65), which holds subject to weak regularity conditions. Observe that applying the K, g 
operator to the function s \-> K, g (h)(s,v) for fixed v, yields 

Kg(K g (h)(;v))(t,s) = K g (h)(s,v)g(t,s) + [ (K g {h)(u,v) - JC g (h)(s,v)) g(du,s) . 



26 



Hence, 



t rt 




K.g(h)(u, v) g(du, s) dB(v) 



K g (Y)(t,s) = / K g (K g (h)(;v))(t,s)dB(v) 
Jo 

rt _ 

K.g, h (s,v,t)dB(v) , 



where we have introduced the definition 



£ g ,/i(s, v, t) = I {v < s} JC g (lCg(h)(-,v))(t, s) + I{„ >s} / JC g (h)(u, v) g(du, s) . 



(59) 



By properties of the Malliavin derivative, see iNualartl (|2006a . p. 25), we have for s < t 



D s JC g (Y)(t, s) = D s I jc g , h {s, v, t) dB(v) = K 9th 
Jo 



(s,s,t) . 



Hence, we find that (after using the stochastic Fubini theorem, see e.g. lNualart &: Zakai dl988L 
Proposition 2.6)) 



Y{s)dX(s) 




t rt _ rt _ 

JCg : h(s,v,i)dB(v)5B(s) + / K. g>h (s,s,t)ds 
10 Jo Jo 

rt r s rt rt 

= / / Kg, h (8,v,t)dB(v)dB( a )+ / / K g>h {s,v,t)dB(v)5B{s) 

Jo Jo Jo Js 

+ / Kg jh (s,s,t)ds 
Jo 

rt rs ^ rt rv 

= / / lCg, h {s,v,t)dB(v)dB(s)+ / / K gA (s,v,t)dB(s)dB(v) 
Jo Jo Jo Jo 

+ / Kg jh (s,s,t)ds 
Jo 

= j j (fcg )h (s,v,t) + Zg th (v,s,tf) dB(v)dB(s) + j JC 9th (s, s,t) ds . 

Hence, we have the explicit representation of an integral of a Volterra process Y(t) with 
respect to another Volterra process X(t) in terms of a sum of a second and zeroth order 
chaos. 

From the general result above we deduce a representation result for J * X(s)dX(s) next. 
Observe that for h = 1 we have Y(t) = X(t). Furthermore, K, g (l)(u,v) = g(u,v). Thus, as 
long as g(u, v) is in the domain of lC g , we find 



Kg{K,g{l){; v))(t, s) = fC g (g(-, v))(t, s) = g(s, v)g(t, s) + / (g(u, v) - g(s, v)) g(du, s) , 



and 



/ IC g (l)(u,v)g(du,s) = g(u,v)g(du,. 

J V J V 



and hence we have 

K 9l i(s, v, t) = 1{ V < S } (g(s, v)g(t, s) + (g(u, v) - g(s, v)) g(du, s) \ + I{ v>s } J g{u, v) g(du, s) 
= ^{v<s} (g(s,v)g(s,s) + / g(u, v) g(du, s) ) + I {v>s} / g(u,v) g(du,s), 
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and K.g t i(s, s, t) = g 2 (s, s) + f g(u, s)g(du, s) . Therefore, we obtain 
tC 9t i(s,v,t) +K. g)1 (v,s,t) 

= ^{v<s} (g(s,v)g(s,s) + J g(u,v)g(du,s) S j + I {v>s} J g(u,v) g(du,s) 

+ 1 { s <v} (g(v, s)g(v, v) + J g(u,s)g(du,v)j + I {s> „ } J g(u,s) g(du,v). 
Altogether, we get 

J X(s)dX(s)= J J (jCg A {s,v,t)+jCg A (v,s,t)^ dB(v)dB(s)+ J K g ,i(s,s,t) ds. 

Let us compare the previous result with ^X 2 (t). By again using integration- by-parts in 
Proposition [21 we have 



* 2(t> ~ (i 



* g(t, v) dB(v)^j = ^ g{t, v) dB(v)^j Qf g(t, s) dB(s] 
g(t,v) ^ g(t,8)dB(a)J SB(v) + J g(t,v)D v (J g(t,s)dB(sfj dv 

rV ft rt rt 

g(t,v) / g(t,s)dB(s)dB(v)+ / g(t,v) / g(t, s) dB(s) 5B(v) + / g 2 (t,v)dv 



/ g(t,v)g(t,s)dB(s)dB(v)+ / g(t,v)g(t,s)dB(s) dB(v) + g 2 (t,v)dv. 
o Jo Jo Jo Jo 

Here, we have used the Fubini theorem for double stochastic Brownian integrals once again. 
Hence, we have 

\x 2 (t)= f [ g(t,s)g(t,v)dB(v)dB(s) + l- [ g 2 (t,s)ds. 
z Jo Jo z Jo 

6 Conclusion and outlook 

In this paper we have developed a stochastic integration theory with respect to volatility 
modulated Levy-driven Volterra (VM.CV ) processes. Our results extend previous work in 
the literature to allow for stochastic volatility and pure jump processes in the integrator. So 
far, we have treated the case of a VMCV process driven by a Brownian motion and a pure 
jump process separately. In future research, we plan to address the problem of accounting for 
a Brownian motion and a jump process simultaneously. 

Moreover, the results in this paper can be regarded as the key building block for developing 
a general integration theory for VMCV processes defined on the entire real line, i.e. for 
processes of type 

X{t)=[ g(t,s)a(s)dL(s), t G R. (60) 



Note t hat stochastic integrals of type (|60p can be defined using the recent theory by IB asse- O'Connor et al. 



(|2012l ) , who studied stochastic integration with respect to so-called increment semimartingales 
on the real line. 
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